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ABSTRACT 


In the heterogeneous flow of finely divided particles and gas 
in a rocket nozzle, the b cb enapetens that the flow enters the nozzle 
with no lags, that the radiative transport of energy between particles 
is negligible, and that there are no temperature gradients within the 

' particles, are commonly made. In this paper, each of these as- 
sumptions is investigated in detail, and they are shown to be reagon- 


able for typical rocket nozzles. 
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I. INTRODUCTION 


In recent years, it has become common to add metals to the 
propellant in rocket engines to increase the chamber temperature. 
Unfortunately, most of these additives form some sort of a solid 
particle in the exhaust,and consequently can cause a depreciation in 
the thrust and specific impulse from the values based on gaseous 
products. That is, instead of an expected 20 per cent increase in im- 
pulse, say, an increase of 12 per cent is noted. 

A very common mixture is a propellant of ammonium per- 
chlorate oxidizer and a polysulfide fuel with aluminum added. For 
the numerical work in this thesis, this propellant is assumed with 
pure aluminum added to it. 

The amount of metal that can be added is limited, since there 
is a requirement for the fuel to act as a binder for the oxidizer and 
metal. Usually, the designer will add as much metal as he can; thus 
a method is needed whereby he can correlate test engine performance 
with the predicted performance. 

Such a method is presented by Rannie in reference 1. In ref- 
erence 1, Rannie used a one-dimensional hydraulic approximation, 
which is known to be valid -_ homogeneous flow, for the heterogeneous 
flow of gas and solid particles. In addition, he assumed ({) that the 
flow entered the convergent section of the nozzle with no lag between 
the gas and solid particles; (ii) that the radiative transport of energy 
between butcher was negligible; and (iii) that the thermal conductivity 
of the solid particles was sufficiently high compared to that of the gas 


that there were no temperature gradients within the particles. 
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Assumption (i) is investigated in Part I] by solving the equae- 
tions of motion in a cylindrical combustion chamber and determining 
the velocity lags at the entrance to the nozzle. Axially symmetric 
flow is assumed in the chamber, and the mixture is approximated by 
a heterogeneous mixture of gas and solid spherical particles. A 
constant source distribution of gas and particles is assiimed along 
the cylindrical wall, and viscous effects in the main flow are ne- 
glected. Turbulence and temperature gradients in the chamber are 
also neglected. It is shown that the particle lag at the entrance to the 
nozzle is indeed very small compared to the gas velocity, and thus 
the assumption that the lag is zero is a valid one. 

The magnitude of the effect of radiative transport of energy 
between particles is estimated by adding a radiation term to the equa- 
tion for the heat balance of a particle and investigating the size of 
this additional term. The ratio of the absorption cross section to the 
geometric cross section is a necessary parameter in this term, and 
a value is not available for aluminum osdde particles of these small 
sizes and high temperatures. However, for any reasonable value of 
this parameter, the radiative transport of energy is shown to be neg- 
ligible compared to the heat transfer by convection and conduction. 

The temperature distribution within a spherical particle is 
found for unsteady heat flux over the surface of the particle in terms 
of the Laplace transform. By expansion of the transform, the solu- 
tion for uniform temperature within the particle is recovered and the 
magnitude of the correction term is determined. The correction term 


is shown to be negligible for practical cases. 
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Ij. TWO-DIMENSIONAL ANALYSIS OF THE FLOW 


IN A COMBUSTION CHAMBER 


Following the notation of reference 1, let c< be the mass 
fraction of particles in the heterogeneous mixture, Py the gas den- 
sity, and f; the density of the solid particles. Then the density of 
the mixture is given by 

| -o& o 


{ 
—_—— = + 
. ¢ = (1) 


Let A be the cross-sectional area of the cylindrical combus- 








tion chamber and m the constant mass flow rate of the mixture. Ap- 
plying the condition that the net flow rate into any volume must be 


zero, the continuity equationsfor the solids and gases separately are: 


2 ](j—«) Pg] . H(t) PY r}=o (2) 


Bel 4s] + 5S See] ee 
‘ (3) 


where Ug and U; are the axial velocity components of the gas and 
solids respectively, and Vg and Vp are the radial velocity come- 
ponents of the gas and solids. 

By considering the momentum transport through a small ele- 
ment and neglecting friction forces, the momentum equation in the 
axial direction can be written 
2. ( \—«)P Ug | Nig & ( -«)PUs V9 r| +2 fap Us | r 

r OX 
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(4) 
os Vert Jeers 
= alia | ax” 
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and the radial momentum equation is 


K(x Pusvslre ofa PV. r| + Ju“) 4g 45 | ie 
plover] arg +P 


Approximating the particles by spheres of radius a , and 


(5) 


assuming that a modified Stokes drag law holds with a correction fac- 
tor {, to allow for deviations from Stokes' flow, the force balance 


in the axial direction on a particle is 





Svat, [US + Ve SMe] = Gta (Ug-l)E — Bora? ger 6) 


Ox fi 2 


The steady state drag law is assumed valid. 4 is the viscosity of 
the gas and is constant since the temperature is assumed constant, 
and the force due to the pressure gradient is added. The cocfficient 
f; is a function of the Reynolds number and the Mach number of the 
gas relative to the particle, and is equal to 1.0 in the Stokes flow re- 
gime. 


A similar equation for the radial force on a particle is 


4 ro3e} U,oVe 4 V; OM] = -V)t 4 gp aidp 
4 ral] Uss¥e + Ve Si] bttaa (Wyle) — 2 agp (7) 


Equations 6 and 7 are expressed more conveniently in the fol- 





lowing form: 
aire 3p 
Ug— Us = 4 7 ¢ {Usee + Us + st} (Sa) 
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and 
Vee a ae a fe $f, (u av ov + $#1 8b 
¢ qu oe ee _ 


Letting U.= Ug + (Us - Ug) and = Va +(Vs-V, in equation 


3 and subtracting it from equation 2 
Safar] flee teyps larly Veo 


Introducing the following non-dimensional variables p= a : 
ies vA \< a4 _ y= a , where VV. is the radial ve 
locity component of the gas at y= | just outside of the boundary lay- 
er and b is the radius of the burning surface. To simplify the 


equations, also let 
Ee = a af Vs Fy 
1M 


Introducing these in the above equations, we obtain for the combined 





continuity equation 


Salta eared Sy Y-4lF W]= ea 


and for the momentum equations 


Hara} efi~ «P+, 4 7] + + [xed] i 
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and for the force balance equations, 











ya OF, ( Jp 
$4 oo = ef &55 ap ae - +t] (13) 
2 Paes ¥ p 
Yg Ys = ¢€ { ?. je + % ots ae 
: eo a4) 


With the help of equations 2 and 3, equations 11 and 12 can be 


further simplified by rearranging so the dependent variables are qo , 


ica aur » and ts— ta 


Thus equations 1i and 12 become 


ec. * ONG ont xP{ Beth + ,— Py SP +{f, yh) 
12. (15) 
+} 14 ol + (¥s 1) ot = (Ys — ¥en" tl = v ‘Ne 


Ye, ol 
Pay he PSE a oy TEs (6-4) 28 Het —. 





Y< ~ oY se) = cae a ) 
+ te er 4) +(s “hea + ['ts-¥y) 57 J ve oh 


(16) 


First Approximation, Zero Velocity Lags. 


When the velocity of the particles approaches that of the gas, 


(i-e., as a. aa O ys ce —> OC 


and >, —> te and > ¥g 





Then 
p= ty =o, and {fs= ‘tg= Pe (17) 
where the subscript zero indicates no lag. 


Then equation 10 reduces to 


OP. O Ye tae 
ee + 
o¥ 24 ” 


and equations 15 and 16 reduce to 





oe (18) 


fo [po Sfe + ¥ Ste] = - 1 se (19) 





O E oO 4 Vee O 
MEG P) Yo = = ale On 
fo | P oS aa O Y | — ven x 


(20) 


It is a simple matter to verify that the solution to these equa- 





tions, satisfying the boundary conditions (a0). = oO at E=O0 
and Y= -—/. at ea | , is 
, Po = rs a 
Y= -7% (21) 
ae + Const 
s= 26 “27 i 
Ve 


Second Approximation, Small Velocity Lags. 


We shall assume for the solution that the variable may be 


expanded in powers of € ; thus 





32 


= P . + € $9, Pa eee 
?. Mis aye co € Ps, pores 
Tq = ie 1 16 Ye, piace 


(22) 
Yo a Yo te Ys toe 
a Sots + 6 Cape: 
a e. + 6 Cee 
The criterion for these expansions to be valid is that 
ga Pe = ae Fe aaa (23) 


and similar inequalities for the other variables. Here € is merely 


an indicator of a smail quantity. Substituting these expansions into 


equations 13 and 14, and equating coefficients of ©€ 


bet, = $28 + ¥ Ste — SHE 





oe oy is 
and 
Pal Dies, ae ie 
ty Ye, = P05 t HOSE 


The density of the mixture, (. 


‘e ; of the solid material, so those terms can be neglected that con- 


tain fo/ ‘ ; therefore 


» ig very much less than the density 


Z O Po 
Pa, _ Po sae Po ote T Yo oT (24) 


wit 


and similarly 
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Phe 
1g, — Ts, = Fo a eis oF (25) 


Now defining the first two terms of equation 10 as F , then 
a Aer 


Equation ¢ can be written as 


_ s[<Pe] - 9 Sl<P x] = 


Then differentiating this by parts, 


= | Sf] “4 Sle tll ~ PS Py, Se + 


or 
a es Ae = 27 
F Cpa) P 4 Se aa =o ‘<a 
Also comparing equations 26 and 10 we see that 
al LOCH 4) | 
ete: | AP LR | 28 
F = 3¢ P($4- #.)| + 7m P(%-%) (28) 
From the first approximation, we may evaluate equations 24 
and 25 as 


pa. = 5, =“ 4% 


(29) 
Gn ts, i a 











=102 

Substituting these into equation 28, the solution of F is 
F= 64,F, (30) 

where use was made of equation 22 and only carried to order epsilon. 


Putting equation 30 back into equation 27 and incorporating 


equation 22, and rearranging, 


25 Se -7 5 = 6 4,(1—%e) 


(31) 
Particular solutions for equation 31] are 
Kd, = 3n,(1-a,.)4n€ +4 (32a) 
kL, = wiki ys 222 ee 7. B 
(32b) 


The boundary condition for “, is chosen so the mass flow 


rates of gas and solids are the same as with no lag. . Since 
4. (=e, Vol le 5) 
° A (82) 
AAs ote fVo} I+ e Ge thr 6.) 
h 
the conditions that 9/4 and "s/4 are the same as with no 
lag, are 
19, =e 
>i + %, we es 
ao fo 


From equation 1, upon substitution of equation 22, and dropping 


= 
terms of order € , 








~lie 


os FoK (35) 
Pesan 





Substituting equations 35 and 29 into equation 34, we find that 


d= —dX,(/—x.) at y= (36) 


satisfies the condition that the last term of equation 34 vanish. 
Clearly the solution 32a is not applicable, and applying condition 36 


to equation 32b, 
=-b to (/-de)/ Lay tt] (37) 


The boundary layer and friction forces along the wall are ne- 
glected and the mass fraction of particles to the mixture is assumed 
to be constant along the length of the chamber since the propellant is 
presumably uniform in mixture. 

This solution does become infinite at the centerline, so the 
solution is valid for the whole chamber except on a line down the cen- 
ter of the chamber. This singularity does not cause any difficulty 
because the mass flow rate in the axial direction is finite, as will be 
shown below. 

Now substitute equation 22 into equation 15, and put f he téP, 
then equate the coefficients of € , dropping known zero derivatives 


(i.e., 240/y and ote /de )3 the result is 








& Pe oR o> 
C #55 5a 41.5 Be + Fotho 5 ale & Y% D Gt 


_ moe $4,-%s,) | eto) yee oP. 
a a y= — Lofel $4,- Fs, OF Vis OF 








and upon substituting equations 21], 29, and 35, 


a © +" 2 $4 +2€ ta ees 15H mee 5 Lok = 1 of, 


| amg favo OF 


Combining the second and third terms and putting in the value (equa- 


tion 37) of %, , we obtain 


24k. (Lng 4 Lead (€43:)- ioe bx. Fo 57 t <p (38) 


Similarly, equation 16 becomes 


= b o9( Ln 7 +4) — aa 


: 
5] +25 +2 KoI= EET (39) 


-2 (y 
3" 


Let us assume a solution of equation 39 to be 
g,.= ~SHen fre (40) 


Substituting equation 35 into the continuity equation for the gas (equa- 


tion 2) leads to 


oFfe, _. LL 2(% 4) a 
UT 
Upon substitution of V4, into equation 41 and integration, 
ty, = bxo Ee Lay tale & (42) 


where the constant of integration disappears since Fa, must be 
zero at ss © ~ 


From equation 29 we find 
1 = EXo¥ Lny eee == tae 
and (43) 


Ts, = ~3el0y-ay -% 








alae 
To verify the solutions, substitute the values for tg a f qi» 
ts 


tions (equations 2 and 3). The mass flow rate of gas for a cylindrical 


eee and f Sf (equations 40, 42, and 43) into the continuity equa- 


chamber of radius b » with no lags, is 


6 


me = ib (1K) Blo te (44) 
and with velocity lags, is 


hg, = ETE (Vas vote | (Ir bt) 1 dy (45) 


6 


Substituting from equation 42 and integrating, equation 45 becomes 
Mag = Ma | (46) 
The mass flow rate of the particles without a velocity lag is 
Pee Th, dol. >, (47) 
and with a velocity lag, 


mM 
=o 


ey } “toe dof [rve(S 7 B+ fa)] yy (48) 


Evaluating the integral 


ay; mee \s (49) 


showing that the mass flow rate of the particles does not change when 
velocity lags are introduced. Hence, the assumed solution and con- 
stants are consistent with mass conservation. 

To obtain some numerical results, a specific chamber must 


be chosen. The following values are consistent with the numerical 








ee 


example given in reference 1. 


—4 
QmeO, 266 fa A= 2.5 * To “aa 
ai 7 | 
Mm = 1.5% 10 lb-Sec/L, wy =0.5 'Ysec 
fe = 210 lbp 3 2 = loco !b/\n* 


1 = Lococo °P X= a, or 


The chamber length is chosen as twice the diameter, and Stokes 


flow will be assumed (f= Io) . Then 


ves Mika Bs eZ 
E = Big nto ie el ae 


The magnitude of the velocity lags in relation to the no-lag 


flow velocity is obtained from equations 29 and 2] 


£9-¢ 26 (50) 


which verifies that the lags are indeed small. 

The uniformity of the flow across the exit plane was assumed 
in the analysis in reference 1], and Table 1] is presented as a verifica- 
tion of this assumption. Values of © i a Fa S go » and 
EI / [ae were calculated as functions of i » the non-dirre nsional 
radius. 

The Reynolds number of the gas relative to the particles was 
computed at the exit plane using the particle diameter as the charac- 
teristic length and found to be 6.5 x 0" ~. The Mach number was 
also found to be 1.36 x lon. Combining these two values, a value 


for M/Re of 2.09 x 10°? is found, which indicates that the particles 


are in the Stokes flow regime, and thus the correction factor fy is 














a5. 


TABLE 1}. 
Y é Pa Peeiat ~ AL. 
Po Fo te 
ak -7.35 x1074 -3.69 x io"? 3,63 x 10% 
ip ~6.10 x 1074 ~2,43 x 107% 1.73 x 107° 
% ~4,75 x 1074 256 x lo h. 26% Yor 
4 “Ade x07" 2051 Six gore 9.25 x 107> 
nc ~4,07 x 1074 -4,21 x 107° 6.49 x 107° 
6 a20ee x 10°~ ~2,29 x 107° 4,24 x 107° 
7 -9. 35% 107" 43.14x 107° 2035 x 100m 
8 ~3.05 x 1074 . Eeno 6,95 x 107m 
9 ~2.81 x 1074 8.41 x 107° -7. 58 x 1074 
1.0 oeese 10-9 1,10 x 1074 ~2.06 x 107* 





approximately equal to one, as taken above. Reference 1 gives a dis- 
cussion on the calculation of the parameter i>: Using equation 66 


and equation 67 of reference 1, 
f= lo+ 459 3- = 1.0096 (51) 
d Ne 


Thus, from the above analysis, the assumption made in refer- 
ence l, i.e., that at the entrance to the nozzle the flow was uniform 


and without lags, is a reasonable assumption. 











-i6- 
Wil. RADIATION BETWEEN PARTICLES 


For a one-dimensional steady flow of a homogeneous fluid 
without heat addition, the energy equation, in terms of internal energy 


per unit mass, is 


PAL (e +) + Api = Const, (52) 


and for the heterogeneous mixture of gas and particles is 


S 


(—«) FAL a (Gyr Wt PAU, (e_r Us) Ap(Ug jal pus Som const. (53) 


where & is the mass fraction of particles (subscript 5 ). This 
equation is general and independent of force interactions and heat 
transfer processes between the particles and the gas. 

If radiative transfer of heat between particles is allowed, an 
additional term must be added on the left hand side of the equation, 
since the equation does not take care of interactions and heat transfer 
processes between particles. The energy transfer between particles 


due to radiation is 





1 SS ee (54) 
i a A 
where 
3 
a - Cea (55) 


may be regarded as an effective conductivity due to radiant energy 
transfer. The coefficient Og is the Stefan-Boltzmann constant, 
and £, is the Roeseland mean free path. Then —— A 
must be added to the left hand side of equation 53 to account for radi- 


ative transfer. 











a7 - 

Equations 54 and 55 are applicable only when Xp is very 
small compared to the length scale appropriate to the problem, for 
instance the nozzle throat radius, at all important wavelengths of 
radiation. 

The variable mass fraction «< in equation 53 can be elimi- 


nated with 





Lo 


obtained from equations 2 and 3. In cases of interest, the volume oc- 
cupied by the particles is negligible when compared to the volume of 


the gas,and the assumption 
q = Cl -X) P (57) 


is valid for a particle material density of 210 Yb /ft” as we have chosen. 
Since the gas density will seldom be greater than 0. 4 ie, this ap- 
proximation is valid over a wide range, but should be checked for any 
particular problem. 

Substituting equations 56 and 57 into equation 53, the energy 


equation becomes 
(1 —4,) €q tos + rq -Ld, o( Ug u 5] ae? 


; (58) 








XoP dTs  _. 
je gal ~ Ke) a 


where h, is the enthalpy in the chamber. 


For a perfect gas with constant specific heats 


ey t =Cylg th Tq =Col (59) 
pet 3 











Ss 
and we now introduce a specific heat, C Sie for the mixture with a 


particle mass fraction ™, , 


Co, = (lao) Cp + Xo € (60) 


and modify the gas constant to 
(Seas (61) 


The right hand side of equation 58 is Cees where a is 
the stagnation temperature of the mixture. Then equation 58 may be 


rewritten as 
{5 a 
(Ke) Cp Tg +ao0Te + £ (1 xe uly + Lao Ue “ai, i-G)ets = Croke eet 


where equations 56 and 59 were used. 
a = 2 a z 


and combining the terms of equation 62, 


Cpe la t 4. +xoC(Ts-Tg)+ 2 a iy" pa = Cpolg (63) 


Neglecting radiation, the heat balance for a sphere of radius 


A and temperature ie is 
4 nas dlkgess ee. 
g WO fe Us T= EE (Ts-F) ae sa) 


where ka is the conductivity of the gas surrounding the sphere and 
Fy is a correction factor that takes into account the influence of cone 
vection. For the Stokes regime of flow, san =1.60 , and it decreases 


at higher Reynolds numbers. 











es 


Rearranging equation 64, 


71 - Us Te (65) 


To express equation 65 in dimensionless form, the following 


notation from reference 1 is introduced 


“L= VR 
im io T (66) 
m= Yr § 


giving 
3g": eee fi 
Te = Teves ioean " % Te (67) 
where A 
ee 
Cee: os P. (63) 


has a numerical value not much different from unity. If the particle 
radius is extremely small, € approaches zero and since the other 
factors on the right hand side of equation 67 are of order unity, 

7, =o and the temperature 136 is zero. The velocity lag is 
also proportional to © so it becomes zero as well. As 2 increases, 
€ becomes larger and the principle of the perturbation procedure in 


reference 1 was to expand Tq ; i: i Pq ,» and Ps in power 


series in € 
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With radiative transfer, there is an additional term in the heat 
balance for the particle. Suppose a particular particle is located at 
X=X, . The radiative transport of energy in the X - direction per 
unit area normal to X inthe space surrounding the particle and close 
to the particle is 


- [k, fi] . fa (k, val X~ Xo) 


o 


Hence the net radiative transport of energy to the sphere is 


[ara sin 0d Cos of (ke -}¢. (a BY (a cose)| 


Ko 
Xo 
and the heat balance for the particle becomes 
ey d io 
3 ra © CUE = — B55 rar Arpad (bh ab. (69) 


which replaces eguation 64. Rearranging and introducing the dimen- 


sionless variables from equation 66 


 ¢ ots dt 


ee ee 
Tg T= Leta dé TBP as a soit de (k) (70) 


The radiative mean free path Lp is given by 


Laz 3 (71) 
x ntar @ 














Zle- 
where n_ is the number of particles per unit volume and @ is the 
ratio ci absorption cross section to geometrical cross section. ro! 
the identity 
> oe 
ngtrals = 7 fy 
and the equation of state for the gas 


afte es. 2a. Som et Ce a 
tine RU Qe 7G 3 7 


Substituting into equation 71 above 


| — ea 
3) Fra Q 
and hence ee from equation 55 becomes 
Y ae ; 
+P es 
c 


Substituting for £, into equation 79 


oA. 


am 
or, 
ty* 


eee RT. a Vy 
eet = ef, 6 3 are Me meats 5 pal $A) 


a 3 2 ColkoTe 


where @ is assumed to have small variation with 3 and can be 
taken out of the differential. Since Ce/e = Y or » equation 76 


can be written as 


Vp ; 
Laer = ef). Bf at ras P. mae ( Pay 77) 


where o- Tk 
















“220 


Since the factors multiplying f, in equation 77 can be con- 





sidered to be of order unity, as is the first term in the square brackets, 
a value of NS of order unity implies that radiative transport is as ef- 
fective as the heat transfer by convection and conduction. In the exam- 
=% 

ple of reference l, { < 3 x |o/@ ; hence, L, is extremely small 
evenif Q is quite emall. Hence the influence of radiative transport 
on particle lags is neplizible. 

Tne zero lag solution is also affectec by radiative transport of 
energy. The energy equation 65 becomes 


Le a ke Croke) Ear 
1 ar) Yo Po aa  GlRele PoC. dé 


Upon substitution of equation 55, we find that the condition for radiation 
to have negligible effect on the zero lag solution is 
4 
Tele Lg a 
cle 
since the non-dimensional terms are all approximately of order one. 
The concept of an effective conductivity has meaning only if 
h 
m , RT —— — 
ER... << | . Sor chemical rockets ““ PMT | = |0 
tlence the radiative term gives negligible contribution. The ratio 
fe/ rr is 
Wigs. (2. Aie Kilegeeeaee 
“~ -S = 
In the example Yr & ZK 10 , (SRTe az SAIC , ame 


i& |,57 . Evenif @ is small, LR po one 


(] 
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a. 
IV. EFFECT OF PARTICLE CONDUCTIVITY 


The effect of finite particle conductivity is investigated by 
solving the heat conduction equation for the particle with an appropri- 
ate boundary condition, then expanding that solution to investigate the 
size of the correction term. 

The heat conduction equation for a spherical, uniform, solid 


particle with spherically symmetric temperature distribution is 


pies @ V ra eee (ret) 


OT eo Gag (79) 
where QO is the diffusivity of the particle material. 
The initial condition is 
T= 19 (eee E—c (80) 


and the boundary condition at the particle surface is 


) O) 
ks (SE) = kg i) es (81) 


O ™~ /v=A-o 


where de is the chamber temperature and C is the time measured 
from the instant the particles enter the convergent portion of the noz- 
zle. The conductivities of the solid and gas are {&, and Kg ree 
pectively. | 

For steady state, with spherical symmetry, the equation for 
heat conduction in the gas is 


(6 5) =e 


ac De (82) 








@Z24e 


which has the solution 


T= 28 
ia 


At r= co , T= Te , andthus B= TJq andat r=a , 
4 4 


T= Ta which gives A=a(Ta-Tq}. Thus 


= 
= = (Tel) 


er ee 


2a 


i. 


or 


kg(SE) 0 = te (eT) 


(83) 
For quasi-steady heat flow in particle and gas, substitute 
from equation 33 into 81 to obtain 
K 5c) bg 
=e m. ~(Tqa-lTa (84) 
Y /psa—o — ( J ) 


This will be used for the approximate boundary condition for the heat 
conduction equation of the particle. The solution to these equations 
(49 ,89 , 84 ) will be obtained by making use of tae Laplace trans- 
form. 


Define 
ple rode =Lf{FOP= FH (85) 


and 
yl = & (36) 
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Then equation 79 becomes 
i 


-£fe+ FP ° R= 
2 1S 2 ao (87) 
e ov —- 
Since there is no source at ~=OoO , S-— oC 8 r=c and 
ete ==, 7 = =. ae 
therefore ] = © where | = oVar=- pete =e 


which means that 


/ 
w~O-6—->o as Yo (88) 


The solution of equation 87 is 


{ 


S6=-Ae + Be + Vor (89) 


Substituting this into equation 88 


i) _ 
as. 3 = +} VE (Ae an = > aa T. 


fe NR pe 
[se acre om 
and to satisfy condition 83 , b= —A » thus 


“Wee ® - ven 2% Sink gr +Ta 
es (91) 


Substituting this into the boundary condition (equation 64 ), 


we obtain for the left hand side 


kK St) = aa 2 Al $2 4r— LY a} 





=OW 











ao 


and for the right hand side 


ff Aa(ty Taf = BLT) AMS Te] 


since | = Tao at Y =a . Combining these equations, equa- 


tion 84 becomes 


ee ea a 


and solving for A 


Aa akg(Ty-Te) a 
Ks] ga cosh ga - Sinh ga (1s) 


Substituting equation 93 into 91 gives 


T= &k Ty-P)| Stine i pT, (94) 
Y Ke 4a coshga—Sink 4 a (1 Sap 


Let us now define an average sphere temperature ile as 


Ee L 
T= i 4frr [dr 
; Beige 

JS 


and hence 


— 2% y= 
Te = Sr| rT de (95) 








yg 


Then ae will be compared with the uniform sphere temperature in 


the heat balance equation for k. infinite (reference 1). 


Cera Cues = ee 4 oh, 
a a le * Ue 7) a (96) 
Aliso define: 


ks Nema Sink ga.(l = “K,) 


so that equation 94 may be written as 


= SEL 2 a Te (97) 


Then substituting and evaluating equation 95 


ap = raat coshga =e ae gs} ee lia (98) 
Bo B= 


Because large values of time TC correspond to small values 


of g , expand the expression for J; ina power series in G to 


obtain the form correct for small Ge » i.e., 


= L hoe ae l : eee ga ms 

,)- 1+ gay + son (ga)? + gal + (Te (99) 
tieyele 
(Erte gaytedh “(4 (5+ Kigatnal ir is Ne 
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Retaining only terms to order (ga)? in equation 99 


75% (Tea ale eehealt eal] 9 soa] 


which reduces to 


na 2a ie k sa : 
U4-Te ai (Ta) TS Steg) - 1 (ga}] (100) 
Recalling that q= = = Bo we find that for 
S 


equal to infinity, equation 100 reduces to 


Tes (4g Seth 


(101) 
and taking the inverse Laplace transformation of this leads to 
Ik Ta~Ts)= us foc nee 
s Gis 
2(T4-T)73 “dt | (102) 


which is identical with the assumption used in the heat balance of a 
particle of infinite thermal conductivity, i.e., equation 96. 

Reference 2 gives values for the thermal conductivity of 
aluminum oxide and its variation with temperature. Measured values 
are given for temperatures up to 1600°%s, and one more point (1800°C} 
is presented from extrapolated data. The thermal conductivity de- 
creases with an increase in temperature up to 1400°C, and then be- 
gins to increase again. The value given for 1400°C (the minimum 


value of k. for this range) is .0125 on . 22 mer 


cm sec "< 








Bt 0 fa 


36,25 >2U, in, 


ft hr — 


To find a reasonable value for the thermal conductivity of the 


gas, the data from the numerical example in reference 1 was used, 


that is, 


Q) 
u 


0.500 BTU/lb Pr , 


6 


nea 15x 10° 1b. 6ec/ft", 


Pr « 0.74 . 

To be most conservative in evaluating ks/ kq » avalue of 
pA was used that would give a high kg . Since “4 varies con- 
siderably with temperature and the Prandtl number is almost constant 
over a large temperature range, the thermal conductivity varies with 
the temperature in the same manneras /(< . Over the temperature 


range of interest, a power law of the following form is satisfactory. 


ab 
ka = Co me, (TG (103) 
rn ae é| 


The maximum k4 would be obtained when less To , thus 


Ke nn tx = 1412 Bre ia 
fT ae fe 


Thus the value of ‘Vk in equation 199 wilt always be 
greater than 25.7, and thus one may neglect 1/5 in comparison to 
ee 

With this simplification in equation 199 and substituting 


SC py for q” ; 
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Ze a a. 
= joe pee eee cece _! fe 4 | 


Taking the inverse of the Laplace transform yields 


ka (TT \ya NU Ys oteets po cad Ts : 
(Ty Ts | 3 fCa C8 ae — (105) 


Without the second term on the right, this cquation is identical with 
equation 101, and the last term is a correction term. To investigate 


the size of the last term, rewrite equation 102 as 


elle Seka. cere (106) 
ie Roar | 4 2 | 


and differentiate to obtain 


apa deme ane | 
Te = «(45 - 45) (107) 


where <= ? "YE oa . Then equation!&5 becomes 
Be lc di; 3 le (108) 
ke (aT hice] St -S( fH) [$B -s5) 


The condition that the approximation used in the heat balance 


equation of a particle, equation 96 , is valid is 





*(dTq _d 
a at) fae | ee (109) 


At ae 








“~ 


o3}e 


For a first approximation, one can say that 416 /dT and 
ee 7 - are of the same order, while still unevaluated. The coef- 
ficient of the numerator is 1.621 x 10°* based on the constants evalu- 
ated earlier. Since day _ — aT can be at most the same order 
as d Ts, and will probably be less, we can say that condition 109 
is satisfied. 

Since this condition is satisfied, equation- 96 is a valid ap- 


proximation which would be in error by much less than one per cent. 
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LIST OF SYMBOLS 


cross section area 


specific heat 


SAP 4] +4 S[f % 1 


Mach numober 
Prandtl number 


ratio of absorption cross-sectional area/geometrical 


cross-sectional area 
gas constant 

Reynolds number 
temperature 

axial velocity component 


radial velocity component 


radius of solid particle 

radius of combustion chamber 

internal energy/unit mass 

Be aan factor 

enthalpy 

thermal conductivity 

Rosseland mean free path 

mass flow rate 

number of particles/unit volume 
pressure; variable of Laplace transform 


rate of heat transfer 
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\%e 


radial distance; burning rate 
time 

axial. distance 

mass fraction of particle 

ratio of specific heats 
indicator of a small quantity 
pressure ratio 
non-dimensional radial distance 
angle (radians) 

viscosity 

non-dimensional axial distance 
density 

diffusivity 

Stefan-Boltzmann constant 
temperature ratio 


non-dimensional axial velocity component 


«8 2B QO eee 


non-dimensional radial velocity component 


Subsc ripts 


b burning surface 
C chamber 

d drag 

g gas phase 

bh convection 

L with lag 








solid particle 


throat 
constant volume 
with small lag 























